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Abstract 

An interesting empirical result in the assembly line literature states that slightly unbalanced assembly 
lines (in the format of a bowl - with central stations less loaded than the external ones) present higher 
throughputs than perfectly balanced ones. This effect is known as the bowl phenomenon. In this study, 
we analyze the presence of this phenomenon in assembly lines with heterogeneous workers (most of 
them, workers with disabilities). For this purpose, we modify an existing model for the assembly line 
worker assignment and balancing problem in order to generate configurations exhibiting the desired 
format. These configurations are implemented in a stochastic simulation model and the obtained re- 
sults are analyzed. The findings obtained here confirm the existence of the bowl phenomenon in such 
situations and contributes to deepen our knowledge in this empirical phenomenon. 
Keywords: Assembly lines, simulation, disabled workers, heterogenous workers, bowl phenomenon. 



1 Introduction 



Assembly lines are productive systems specially useful for large scale standardized manufacturing. 
The main rationale behind such systems is labour division. This is done by dividing the tasks to be ex- 
ecuted among a number of workers or workstations (Scholl 1999). In each workstation, a subset of the 
tasks is executed. The most classical assembly line balancing problem is known as SALBP (simple as- 
sembly line balancing problem) and considers, among other hypotheses, that all workers/workstations 
are equally efficient. A large amount of research has been done on the SALBP and includes both clas- 
sical and recent surveys and research articles (Salveson 1955, Tonge 1961, Baybars 1986, Ghosh & 
Gagnon 1989, Scholl 1999, Scholl & Becker 2006, Boysen & Fliedner 2008, Scholl et al. 2010). 

SALBP's hypotheses are rarely valid in practical contexts. This fact has motivated the study of a 
large number of its variants (Boysen et al. 2007, Boysen et al. 2008, ?) . In particular, this research 
addresses the case of assembly lines in sheltered work centers for disabled. In this case, due to workers 
heterogeneity, task execution times depend on the worker they are assigned to. This gives origin to 
a problem in which, besides assigning tasks to workstations, there is also the need to assign workers 
to workstations. This situation has been first modelled by Miralles et al. (2007) and named ALWABP 
(assembly line worker assignment and balancing problem). 

Both SALBP and ALWABP assume deterministic execution times for each task and, hence, the 
line's throughput is dictated by the most loaded workstation. In practice, however, task execution 
times are usually stochastic. This may occur due to small differences among parts, changes in worker's 
behaviour, technical issues and many other random events. The randomness in execution times may 
lead to interesting effects such as the well known bowl phenomenon, in which slightly unbalanced lines 
in the format of a bowl, i.e., with central workstations less loaded), or with central workstations with 
more stable loads (less variable execution times), have a better productivity (Hillier & Boling 1966). 

The theory behind the Bowl phenomenon is not yet fully understood. Nevertheless, the effect may 
lead to adjustments in optimization methods (heuristic and exact) in order to enhance throughput of 
assembly lines in practical contexts. Thereunto, simulation models are helpful, because they allow the 
modelling of execution times with higher precision. 

In this study, an assembly line simulation model for productivity analysis is proposed. In particular, 
the interest is in the ALWABP case, for which no similar study has been developed to the best of the 
authors knowledge. The main goal is to verify the bowl phenomenon in assembly lines with heteroge- 
neous workers. For this purpose, a mixed-integer programming model capable of supplying solutions 
with unbalanced bowl-shaped profiles (with greater or lesser slope) is proposed. These solutions are 
used as parameters in the simulation model, also developed in this work, and the results are discussed. 

The remainder of this article is organized as follows. In the following section, the ALWABP prob- 
lem and model are presented. Also in this section, this model is modified in order to generate solutions 
with unbalanced bowl-shaped profiles. Then, Section[3]focus on the proposed simulation model devel- 
opment while Section [4] presents the computational experiments and results. General conclusions end 
this paper in Section|5] 



2 ALWABP 



Miralles et al. (2007) have proposed the ALWABP motivated by the situation found in sheltered 
work centers for the disabled, in which assembly lines are operated with disabled workers. This prob- 
lem has fomented an intense amount of research both on the original problem (Miralles et al. 2008, 
Chaves et al. 2009, Moreira et al. 2012) and on some of its variants (Costa & Miralles 2009, Cortez & 
Costa 2011, Araujo et al. 2012, ?). 

In order to model the ALWABP, let N be a set of tasks. The execution time of each task i G N 
depends on the worker it is assigned to and is assumed deterministic. Indeed, p w i is the executing time 
of task i by worker w. Tasks i G N must respect a partial ordering. We use the notation i <j to indicate 
that task i must be completed before task j is initiated. Also, let W represent all workers available and 
S represent the set of workstations. Defining binary variables y sw (equal to 1 if, and only if, worker 
w G W is allocated to station s G S) and variables x sw i ( equal to 1 if, and only if, task i G N is 
executed at station s G S by worker w G W) and a variable C representing the cycle time (execution 
time of the most loaded workstation), the ALWABP can be written as: 

(1) 

Vi G N, (2) 

G W, (3) 

Vs G S, (4) 

Vi,jeN,i<j, (5) 

Vs G S, (6) 

VseS,VweW, (7) 

Vs G S,Vw G W, (8) 
Vs G Syw G W,Vi G N. (9) 

The the goal of model ([T])-(|9j> is to obtain worker and task assignments to the workstations so 
that the cicle time is minimized. Constraints (|2]) force every task to be executed, while ([3]> and (|4]) 
guarantee that each worker is assigned to a workstation and each workstation receives a single worker, 
respectively. The task execution partial ordering is respected due to constraints ([5]>, while (|6]> define 
the line's cycle time. Finally, constraints ^ make ensure that a task is executed by a worker in a given 
workstation only if the worker is assigned to the workstation. 
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2.1 The bowl phenomenon and the ALWABP 

Manly in non-automatized assembly lines, where people, not machines, are responsible for doing 
the tasks, the execution times are not deterministic. Thus, production rate depends on the distribution 
dictating the task times. It was realized (Hillier & Boling 1966) that if the operation times (task 
times) were following an exponential distribution, then the assembly line throughput could be enhanced 
by introducing specific small load unbalances. This effect became know as the bowl phenomenon, 
since the graphic of total time of each workstation along the line has the shape of a bowl, with times 
increasing from central station(s). Figure [TJ shows an example of such configuration in a line with 8 
workstations. 
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Figure 1 : Example of a bowl profile. 

There is an equivalent of the bowl phenomenon which happens regarding buffers capacity placed 
between the workstations. An assembly line's throughput can be enhanced by placing buffers with 
a slightly larger capacity between central stations (Hillier 1993). Hillier (1993) presents numerical 
results, expanding the findings from (Hillier 1979), which were only available for small lines. These 
results show that, for many cases, with bowl allocation, there can be nearly a 1% raise in line's through- 
put, up to a 2% enhance. This result gradually decrease with the increase of buffer capacities and with 
the reduction of the variation coefficient C, defined as C = ^ , in which fj, is the task execution time 
distribution mean and a its standard deviation. Finally, Hillier (1996) also analyzed the robustness of 
the effect, concluding that bowl-shaped lined provided better results than complete balanced assembly 
lines even with relatively large deviations with respect to the "ideal" bowl format. 

In this work, we try to conclude the existence of the bowl effect in lines with a heterogeneous 
workforce, in order to do so, the first step is to generate solutions respecting the bowl load profile 
(Figure [TJ). This was done in two ways with two different objectives. This first strategy generates 
bowl-like solutions in order to study the phenomenon considering workstation load. The second one 
also targets the obtention of the bowl effect but focus on the execution times deviations: tasks with less 
variable times are assigned to central stations. These two approaches are presented in the following 
subsections. 

2.1.1 Load unbalance 

Constraints ([6]) are appropriate for the deterministic scenario because they induce a uniform load 
distribution among all workstations, hence minimizing the cicle time C. In order to obtain bowl-like 



solutions, these restrictions are modified as follows: 
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in which a s are parameters whose values respect the following equations: 

a\ > 0, 

a s = a\ s \_ s+1 , s = 1 . . . [|5|/2J. 
a s = {3a a -\, 



s = 2...L|5|/2j -1. 



(10) 



(11) 
(12) 
(13) 



Constraints (111 guarantee positive loads for workstations while constraints ( 12 1 induce symmetry in 
the load configuration. Finally, constraints ( [T3] ) define the bowl slope level which is controlled by 
parameter j3, < j3 < 1. Figure [2] shows examples of resulting profiles for some values of a and (3. 
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Figure 2: Load configurations for different parameters a and j3. 



2.1.2 Deviation unbalance 

Other characteristic that may induce a bowl effect is related to task times deviations. Let c s = £ 
indicate the level of variability in a given station execution time s G S. An inverse bowl profile is 
obtained by maintaining the station loads fixed and reducing the task deviation times in central stations. 



The variation was computed considering solutions without restrictions ( 10 1 and by artificially adjusting 
task deviation times from the base case. 

For five stations, the balanced case (c s = 10, s = 1, .., 5) was modified by increasing parameter C3 
by a unit. This change was conpensated with a 0.5 reduction in each extreme workstation, keeping the 
total load and simmetry of the line unchanged. Similarly, for six stations C3 and C4 have a unity raise 
and c\ and cq have a unit reduction; and, for seven workstations, C4 has a unity raise and c\ and cj have 
a 0.5 reduction. 



3 Simulation Model 

Simulation is a technique for performance and reaction analysis of a system. It can be useful both 
in planning stages or when the system is already in operation if careful validation on practical contexts 
is effected (?). Most benefits of simulation appear in contexts in which the system under study has 
many features that can not be modeled properly with the use of other analytical techniques due, for 
example, to its complexity. 



When building a simulation model, it is important to define the inputs and outputs so the results 
may be analyzed. This analysis depends on the nature of the model. If the model is deterministic, 
for example, the model only have to run once for each input one may want to analyze; on the other 
hand, in a stochastic model, it is possible to use a Monte Carlo method: the essence of it is to use 
stochastic distributions, each one representing a process in the model. Then, a statiscal sampling is 
used to summarize the result for each input. In other words, it can be used to analyze process which 
have a degree of uncertanty (randomness) in the variables (Kroese 2008). 

The simulation model built in order to evaluate assembly lines in this work has the number of 
products which left the line as a state-variable, it is of continuous time, discrete state, stochastic, linear, 
dinamic, open and unstable (Ross 2006). The system has been modeled with the commercial simulation 
general-purpose software Goldsim (Goldsim 2010). A schematic representation of the model build is 
presented in figures [3] and [4j which show the external control (responsible for controlling the flow 
of unities between stations) and the internal control (responsible for defining the time the unity is 
processed in each station), respectively. 
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Figure 3: External control. 
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Figure 4: Internal control. 
The external control ensures that there is always a piece available for the first station, this is done 



with the event new product which draws unities from an infinite buffer (buffer DA). The product is 
put into unit buffers so it can be passed from a previous station to a following one. It is important to 
state that these are not real buffers: the product is put in and removed from them in the same moment 
(the buffers are used only for implementational issues). The final buffer number jof ^products stores the 
number of produced unities during simulation. Elements Machine 1, Machine! and Machine^ represent 
workstations and are there to encapsulate the internal logical control of each station. 

The internal control was constructed in the following way: there is a delay (machine. delay 1) whose 
value is the sum of individual task execution times (which depend on the worker assigned to the sta- 
tion). This delay will determine how long the piece will stay in that station. The workstation has a status 
(machine. ready 1 ), which is either ready, in case it is ready to receive another piece, or busy, if the sta- 
tion is processing a piece. There are also three possible events: pick from buffer (pickjromJbufferJ), 
put in buffer (put An -buffer _i) and delay (machine. delay 1). When a product/piece is withdrawn from a 
previous buffer (pick from buffer event) the status is changed from ready to busy and a delay event is 
generated. After the delay has been processed, if the following workstation has a ready status, an event 
put in buffer is generated, otherwise the event waits until the following station is ready to process it. 

As mentioned, the total value of the delay (machine Aelayl) depends on values of task time events 
associated to it. In this study we use a stochastic simulation and these elements are configured with 
two parameters: mean and standard deviation. Both are obtained from solutions generated with the 
mathematical model presented in Section [2] The solution is schematically depicted in Figure [5] 
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Figure 5: Scheme for obtaining computational results. 

The mathematical model's inputs are: the precedence graph, which determines the order in which 
tasks must be executed; workers, indicating how long each worker takes to accomplish each one of 
the possible tasks; and the number of available workstations. The model's solution is a possible as- 
signment of workers and tasks to workstations. This allocation leads to simulation parameters, more 
specifically, the values (mean and standard deviation) of taksJime elements which are summed in a 
machine Jime element. With these parameters, the simulation model may compute the assembly line 
stochastic throughput. 

4 Computational results 

The computational results were obtained with a series of stochastic simulations using a Monte 
Carlo method. Each input was ran 250 times in order to obtain a statistically relevant results. Execution 



times were assumed to follow a Normal distribution, N(fi, a 2 ), with mean jj, and variance a 2 , ji being 



the deterministic value and a 



-j^j in simulations in which the mean was the variable of interest. 



In order to obtain results related to the mean, parameter j3 from equation ( 13 1 was slowly reduced 
in central stations until a solution became different from the previous one, this characterize a new 



depth bowl level to be analyzed. Parameter a s from equations ( 10 1 were chosen so that a more strict 
time restriction in central workstations implies more flexible times on external workstations. For each 
combination of a s and 0, the modified mathematical model was solved. The resulting solution defined 
tasks and workers assignments in the simulation model. The results, considering line productivity for 
each level studied, are shown in Table [T] These tests were done using classical instances found in the 
literature (Chaves et al. 2009). The best solutions of each case are bold marked. Table [2] discriminates 
the load for each station in the best solutions, shown in Table [T] and the corresponding throughput 
enhancement. 



Table 1: Results for unbalances on mean loads. 





Unbalance 


Instance 


level 


level 1 


level 2 


level 3 


level 4 


level 5 


level 6 


HESKIA 01 


144.310 


144.580 


143.750 


142.870 


142.890 


139.270 


138.190 


HESKIA 1 1 


78.804 


73.572 


78.760 


79.016 


78.864 


78.596 


77.248 


HESKIA 21 


66.004 


66.960 


66.012 


65.696 


65.696 


65.140 


63.916 


HESKIA 31 


64.592 


65.076 


63.968 


65.280 


64.832 


65.096 


64.272 


HESKIA 41 


389.250 


390.040 


390.020 


395.710 


403.460 


392.000 


383.900 


HESKIA 51 


275.010 


271.500 


276.760 


276.880 


270.840 


275.920 


276.100 


HESKIA 61 


203.700 


204.660 


206.460 


204.670 


205.910 


205.520 


203.380 


HESKIA 71 


150.530 


147.840 


149.510 


148.660 


148.660 


149.350 


145.860 


ROSZIEG01 


712.790 


704.490 


687.430 


593.570 


677.650 


675.470 


646.990 


ROSZIEG 1 1 


475.780 


457.830 


420.130 


420.250 


420.120 


420.070 


420.080 


ROSZIEG 2 1 


495.350 


459.630 


459.490 


445.550 


459.590 


431.840 


431.700 


ROSZIEG 31 


456.610 


456.490 


456.260 


456.740 


456.570 


429.640 


419.140 


ROSZIEG 41 


1371.300 


1406.900 


1411.700 


1402.500 


1376.300 


1406.000 


1369.200 


ROSZIEG 51 


1251.300 


1286.000 


1277.000 


1276.200 


1286.000 


1286.900 


1277.500 


ROSZIEG 61 


854.430 


874.370 


832.770 


832.550 


832.800 


826.340 


832.800 


ROSZIEG 71 


914.640 


937.640 


920.960 


935.400 


868.460 


872.100 


880.850 



It can be seen (Table [TJ that balanced solutions (level 0) are only the most efficient configurations 
in a few cases (instance 71 from the Heskia family and instances 01,11 and 21 from the Roszieg fam- 
ily). For the other instances, different levels of unbalance presented the best results. On the other hand, 
Table [2] shows that the obtained solutions are not perfect bowl-shaped, this happens because tasks are 
not divisible so a smooth bowl-shaped solution cannot be guaranteed. Examples of good and bad bowl 
shaped formats can be found in instances Heskia 11 and Heskia 61, respectively. Nevertheless, con- 
firming the robustness of the bowl effect predicted by Hillier (1996), both cases presented productivity 
improvements with respect to the balanced case. 



Table 2: Station loads for optimal results. 
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The best enhancement in throughput was obtained in instances Heskia 41, in which a 3.65% in- 
crease in throughput could be observed in comparison to the standard balanced case. Figure [6] shows 
how throughput (vertical axis) evolves as a function of the level (horizontal axis) for Heskia 41; the 
points were connected to emphasize the behaviour around the best solution obtained. In general, the 
curve smoothness cannot be guaranteed due to the indivisibility of tasks to be allocated. 




Figure 6: Throughput graphic as function of bowl depth for mean. 



The results showing how the assembly line is affected for cases in which central execution times 
are less variable than those in the line's extremes are in Table [3] Considering instances with tasks with 
greater deviation on peripheral workstations, the bowl phenomenon appeared not as often as in the 
previous case. Nevertheless, the effect could be observed in some instances. 



Table 3: Results for unbalances in U values. 





Unbalance 


Instance 


level 


level 1 


level 2 


level 3 


level 4 


level 5 


level 6 


HESKIA 01 


144.310 


144.010 


143.560 


142.810 


141.800 


140.390 


138.140 


HESKIA 1 1 


78.800 


78.570 


78.350 


78.000 


77.500 


76.692 


75.658 


HESKIA 21 


66.000 


65.980 


65.860 


65.550 


65.120 


64.560 


63.708 


HESKIA 31 


64.590 


64.600 


64.470 


64.210 


63.896 


63.428 


62.724 


HESKIA 41 


389.250 


388.450 


387.460 


386.040 


384.120 


381.510 


377.130 


HESKIA 51 


275.010 


274.480 


273.810 


272.860 


271.520 


269.570 


266.500 


HESKIA 61 


203.700 


203.750 


204.280 


204.660 


204.910 


204.580 


203.220 


HESKIA 71 


150.530 


150.410 


150.390 


150.190 


149.630 


149.000 


147.580 


ROSZIEG01 


712.790 


712.560 


712.390 


712.190 


712.020 


711.760 


710.080 


ROSZIEG 1 1 


475.780 


475.800 


475.820 


475.870 


475.790 


475.790 


475.680 


ROSZIEG 21 


495.350 


494.860 


493.820 


492.800 


490.320 


487.000 


481.970 


ROSZIEG 31 


456.510 


455.810 


454.900 


453.620 


451.800 


449.470 


445.750 


ROSZIEG 41 


1371.300 


1373.100 


1374.400 


1375.300 


1375.400 


1374.300 


1370.100 


ROSZIEG 51 


1251.300 


1252.400 


1252.800 


1252.500 


1250.800 


1246.600 


1238.000 


ROSZIEG 61 


854.430 


853.700 


852.200 


849.800 


846.400 


841.170 


832.020 


ROSZIEG 71 


914.640 


916.060 


917.200 


918.100 


918.760 


919.120 


918.880 



Table [4] details the line configuration of the best solution found for each instance (associated with 
the bold value in each line of Table [3]>. The table lists the load of each station and the value of the 
U = fx/a parameter (determined according to section [2]). The greater enhancement in throughput was 
0.6% comparing to balanced case. This happened in a 4 unit unbalance, i.e., with parameters Us and 
U4 (central stations in Heskia 61) raised by four units. 

Heskia 61 

206 
205 
204 
203 
202 
201 
200 
199 
198 
197 
196 

01234567 



Figure 7: Throughput graphic as a function of bowl depth for standard deviation. 

Figure[7]depicts how assembly line's throughput evolves related to level of unbalance for HESKIA 
61. The points were connected to emphasize function behaviour around the best solution found. This 
curve tends to be smoother than the one presented on Figure[6]because raising the unbalance, parameter 




Table 4: Detai 


s for U unba 


ance 




Instance 




SI 


S2 


S3 


S4 


S5 


S6 


S7 


Increase 


HESKIA 01 


Locid 

u 


93 

y *J 

10 


94 

y^ 

10 


92 
10 


94 

y^ 

10 




- 


- 


00 % 


HESKIA 1 1 


T C\ 'Af\ 

u 


1 68 
10 


1 68 
10 


1 69 
10 


1 69 
10 


- 


- 


- 


00% 

\J.\J\J /V 


HESKIA 21 


Loctd 

u 


200 

Z. \J \J 

10 


200 
10 


199 

J. y y 

10 


200 
10 


- 


- 


- 


00% 


HESKIA 31 


T O 'A H 

u 


204 
10 


203 

sL\jy 

10 


204 
10 


204 
10 


- 






00% 


HESKIA 41 


T naH 

L^iXJtXKJ. 

u 


33 
10 


35 
10 


34 
10 


31 

^ J. 

10 


29 
10 


35 
10 


35 
10 


00% 


HESKIA 51 


T oaH 

u 


.J 1 

10 


48 
10 


40 
10 


33 
10 


36 
10 


40 
10 


49 
10 


00% 

\J.\J\J /V 


HESKIA 61 


T oaH 

i_/Uau 

U 


60 
6 


61 
10 


62 
12 


66 

VJVJ 

14 


66 
12 


59 

*j y 

10 


57 
6 


59% 


HESKIA 71 


T O 'A H 

u 


QO 

y\j 

10 


61 
10 


91 

y i 

10 


66 
10 


82 
10 


86 
10 


84 
10 


00% 


ROSZIEG01 


T C\ 'A H 

u 


17 



20 



17 



16 





- 


- 


00% 


ROSZIEG 1 1 


Lo&d 

u 


29 
1 


% 
13 


30 
13 


19 

J. y 

i 


- 


- 


- 


02% 


ROSZIEG 2 1 


T oaH 

u 


25 
10 


27 
10 


28 
10 


28 

Z.O 

10 


- 


- 


- 


00% 

\J.\J\J /V 


ROSZIEG 31 


Loctd 

u 


31 

6 


30 
10 


26 
12 


23 

Z,,J 

14 


- 




- 


00% 


ROSZIEG 41 


T oaH 

u 


7 

1 

6 


10 
10 


10 
14 


q 

y 

14 


10 
10 


q 
6 


- 


30% 


ROSZIEG 51 


Load 
U 


6 
5 


11 
10 


10 
15 


11 
15 


11 

10 


10 

5 




0.11% 


ROSZIEG 61 


Load 
U 


12 
10 


14 
10 


16 
10 


15 
10 


16 
10 


16 




0.00% 


ROSZIEG 71 


Load 
U 


11 

5 


9 
10 


14 
15 


15 
15 


15 
10 


12 
5 




0.50% 



U, does not imply a possible task reallocation; this happens in the mean case because it is more closely 
related with the solution obtained from the mathematical model which is restricted by parameter a of 
constraints (fTOl). 



5 Conclusions 



This study introduced a mixed-integer mathematical model to obtain slightly unbalanced load pro- 
files for assembly lines with heterogeneous workers. These profiles were used in a simulation model to 
verify that the efficiency of the lines could be improved if they were slightly unbalanced. These results 
expand the scope of the well-known bowl phenomenon, which is empirically know since the sixties 
but whose theoretical causes are yet not very well explained. We believe this simple study somehow 
reopen some questions on this classical phenomenon and encourages new research on its effect on lines 
with heterogeneous resources. 
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